Interactions and correlations of particulate inclusions in a columnar phase 
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We calculate the elastic field mediated interaction between macroscopic particles in a columnar 
hexagonal phase. The interaction is found to be long-ranged and non-central, with both attractive 
and repulsive parts. We show how the interaction modifies the particle correlations and the col- 
umn fluctuations. We also calculate the interaction of particles with the topological defects of the 
columnar phase. The particle-defect interaction reduces the mobility of the defects. 
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I. INTRODUCTION 

Macroscopic particles suspended in an isotropic fluid - 
colloidal dispersions - have been the subject of intensive 
study for over a century Recently, a novel class of 
materials have been realised in which colloidal particles 
are dispersed in anisotropic fluids, typically thermotropic 
and lyotropic liquid crystals [|[ ||, |4j. The combina- 
tion of anisotropy and long-range order inherent in these 
mesophases leads to totally new kinds of interactions be- 
tween the dispersed particles, which have no analogues 
in conventional colloidal dispersions || ||, Q, g, §]. 

The most important of these new interactions results 
from the elastic deformation induced by a particle in 
the liquid crystalline host. The overlap of the the elas- 
tic deformations of two particles produces an effective 
long-ranged interaction between them. This mechanism, 
which relies only on the mean distortions of the order 
parameter, implies that other sources of deformation like 
topological defects also interact with particles through 
long-range forces. Fluctuations of the order parameter 
together with the boundary conditions on the particle 
surface lead to another novel kind of Casimir-like inter- 
action [fiof , which can often be important near special 
points in the phase diagram. In particle-doped nematics, 
for instance, director fluctuations lead to a Casimir force 
which is stronger than the van der Waals force near the 
isotropic-to-nematic transition Q. Interparticle interac- 
tions are further complicated by the inevitable presence 
of topological defects surrounding the particles M . These 
defects compensate the mismatch in the boundary condi- 
tions imposed on the order parameter at the particle sur- 
face and at infinity. To clarify with an example, a particle 
in a nematic with director-normal boundary conditions 
on its surface is equivalent to a radial hedgehog with unit 
positive topological "charge", and is incompatible with 
the requirement of an uniform director configuration at 
infinity. The mismatch is relieved by the spontaneous 
appearance of a defect of unit negative "charge" , which 
is tightly bound to the particle. Such companion de- 
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fects usually lead to rapidly decaying, repulsive, power 
law forces between the particles The interactions 
mentioned above, when combined with the multitude of 
interactions already present in conventional colloidal dis- 
persions ][lj, can drive the particles into spatially organ- 
ised structures - chains q, jl3fl, clusters 0, and even 
crystalline arrangements 1 12] have been observed for par- 
ticles in nematic hosts. 

As a result of these novel interactions, the static and 
dynamic properties of (with an obvious extension of ter- 
minology) liquid crystalline colloids are quite unusual. 
The change in the mechanical properties is particularly 
dramatic - particles added in volume fractions as low 
as several percent can totally alter the rheological be- 
haviour of the host phase [|[ [|, When the par- 
ticles have a magnetic moment interesting ferro-liquid 
crystalline phases are formed ][l5| [l6| . 

Numerous theoretical studies have been done on these 
and related effects in nematic colloids ]Tt| . Compara- 
tively less work has been done on particle-doped lamel- 
lar phases (T|[ It is therefore surprising tha ; 
in spite of recent experiments fl5| , |l6| , there has been 
no theoretical study of the properties of particle-doped 
columnar phases. In the following sections of this paper, 
therefore, as a first step towards a more detailed theory, 
we calculate the elastic field mediated interparticle inter- 
actions (section flj) , as well as the interactions of particles 
with point topological defects (section III ) . We conclude 
the paper with a critical discussion of our results and 
possible extensions of the present work. 



II. PARTICLES IN A COLUMNAR 
HEXAGONAL PHASE 

Let us recall that the columnar phase is macroscop- 
ically equivalent to a regular array of tubes, aligned 
in, say, the z direction (directions in the x-y plane are 
denoted as _L) and free to slide along their axes. A 
cross-sectional slice normal to the tubes shows a two- 
dimensional crystal structure, which, for the case we 
study, is hexagonal. Mechanically, the phase is liquid- 
like in the z direction, and solid-like and isotropic in the 
_L directions. Elastic deformations are described by a 
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pair of "broken-symmetry" variables u = (u x ,u y ) cor- 
responding to the local displacements of the tubes from 
their equilibrium position. To lowest order in gradients, 
the harmonic free energy density includes contributions 
from in-plane bulk and shear deformations, and out-of- 
plane column curvatures [p3|: 



forces on their surroundings |L9], |39| . Henceforth, we shall 
assume a to be positive. It is worth emphasising that our 
theory is valid for separations large compared to both d 
and a, which represent, respectively, the scales beyond 
which the elastic and concentration field descriptions are 
meaningful. 
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VyUyf + {VyU x + V X Uyf}. (l) 



Here B and C are first-order elastic constants (with di- 
mensions of energy density) corresponding to bulk and 
shear deformations respectively, while K 3 is a second- 
order elastic constant (with the dimension of force) cor- 
responding to column curvature. An elastic description 
of this sort is valid at length scales large compared to the 
column spacing d. The colloidal particles are modelled as 
monodisperse spheres of radius a and are described by a 
concentration field c(r) = J2v ~ r v)> where v v is the 
center-of-mass coordinate of the vth particle. A descrip- 
tion of this sort is valid at length scales large compared to 
the particle radius a, when in effect, the particles can be 
treated as point masses. In this dilute limit, where the 
typical distance between particles is much larger than 
their size, all short-range interactions can be ignored, 
and the particles need to be assigned only an ideal gas 
entropy. As in the theory of point defects in three dimen- 
sional crystals, it is assumed that the particles behave as 
force dipoles ("centers of dilation or compression") f2~if| 
and produce a body force proportional to Vc(r) at the 
point (xi, z) = r. This is equivalent to a coupling of the 
form c(r)E(u) in the free energy, where E(u) = Vj_ ■ u is 
the local dilation. It is not difficult to see that to linear 
order in u and lowest order in gradients, this is the only 
possible symmetry-allowed coupling of the concentration 
to the column displacement |2f|. With these consider- 
ations, then, the free energy density (in units of &bT) 



T{c, u) = Thex + c(r) 



co 



ac(r)E(u). (2) 



The second term in (Q) accounts for the entropy of the 
particles whose mean concentration is Cq. a is a coupling 
constant, which on dimensional grounds must be of order 
Bet 3 . Note that in our units B has the dimensions of 
inverse volume and that a is dimensionless. The micro- 
scopic details of the boundary conditions on the particle 
surface are incorporated phenomenologically in the sign 
and magnitude of a. For positive a the columns are re- 
pelled away from the particle surface leading to a local 
dilation; this situation is most appropriate for inert col- 
loidal particles, whose size a is appreciably larger than 
the column spacing d. For negative a the columns are 
are attracted towards the particle leading to a local com- 
pression; this situation may be appropriate for small bio- 
logical inclusions, which are known to exert local inward 



Effective Interactions 



Integrating out the displacement field u from F[c, u] = 
/ d 3 r J-(c, u) we obtain an effective free energy 



F r = -In / Vue 



-F\c,u] 



(3) 



for the concentration fluctuations. Since our interest lies 
only in the mean displacements of the columns, the Gaus- 
sian integral in (JsJ) can instead be replaced by the min- 
imum of F[c, u] over u. Setting the variation 6F/Su 
to zero, we obtain a pair of coupled linear differential 
equations for the components (u x ,u v ). Transforming to 
Fourier space (ie multiplying by e _lq r and integrating 
over d 3 r), and solving the resulting algebraic equation 
for the Fourier components, we have 



(B + C)ql + K 3g] 



(4) 



This equation gives the (mean) displacement field for an 
imposed concentration fluctuation c q . Using it to elimi- 
nate the displacement field from (||) , we obtain the effec- 
tive free energy as 

F c = f d 3 r c(r) - ll + \ [ c q V(q)c_ q (5) 

J c Z J 



where 



v(q) 
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(B + C)ql + K 3 qi 



(6) 



and we have used the shorthand J q = J . The sec- 
ond term in (||) represents the effective interaction medi- 
ated by the column displacements. The kernel V(q) can 
be identified with the Fourier transform of the interparti- 
cle potential V(r), with the understanding that the cor- 
rect physics is captured only for lengths large compared 
to the particle size. With that in mind, the interparticle 
potential is given by 

V(r) = ^!viG(x ± ,z) 



a 2 AL„ 



dt 



g— z /4:X^x± cosh t 

(X 3 x± cosht) 1 / 2 



(7) 



where A§ = K 3 /(B + C). The Fourier transform G(xj_, z) 
of the kernel (q\ + A^) -1 can be obtained explicitly in 
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FIG. 1: Plot of the equipotentials of V(r) ~ V±G(r). The 
coordinates are in units of A3. The contour lines have been 
drawn on equal logarithmically spaced intervals from .001 to 
1000, and -0.001 to -1000. The potential is even in z. 



terms of modified Bessel functions p9] , but we prefer the 
integral representation given above. 

The potential V(r) presents several interesting fea- 
tures. It is non-central (not a function of r only) and 
radially non-monotonic. This last property prevents us 
from inferring the sign of the radial component of the 
force from that of the potential. Therefore, the attrac- 
tive or repulsive character of the interaction cannot be 
correlated with the sign of the potential, and the forces 
have to be considered explicitly. However, in the z and _L 
directions the potential is monotonic and decays with a 
power law. Near the z-axis (for z 2 ^> 4:X 3 x±) the poten- 
tial vanishes as |z| -5 , leading to a rapidly damped and 
repulsive force along the columns. Near the x-y plane (for 
z 2 <C 4X 3 x±_) the potential decays more slowly as x^°^ 2 , 
implying a repulsive force in the _L direction. In a general 
direction the nature of the interaction is best visualised 
through the plot of the equipotentials of V(r) shown in 
Fig.|l|. From the plot it is clear that the radial compo- 
nent of the force is positive near the z-axis and the x-y 
plane, and negative in the region bounded by the two 
paraboloids z 2 ~ Xx±. Broadly speaking, the force felt 
by a particle at (x_l, z) due to another particle at the ori- 
gin varies from repulsive, to attractive, to repulsive again 
as the polar angle 9 — tan^^z/xjj varies from to tv/2. 
Since the potential is even in z, the sequence repeats as 
the angle decreases from it to ir/2. A qualitatively simi- 
lar angular dependence is also found in the quadrupolar 
forces between particles in a nematic Q; the distance 
dependence, though, is totally different. 

Given the nature of the interparticle forces described 
above, we might expect an attraction induced aggrega- 
tion of the particles. However, repulsive forces due to 
the companion topological defects mentioned in the intro- 



duction may dominate at short distances, and arrest the 
aggregation. In that case, it is likely that particles will 
spontaneously organise in chains oblique to the columns. 
This process will be dominated by Brownian motion, and 
not the elastic interaction, which though long-ranged, de- 
cays too rapidly to affect the kinetics of aggregation. 

The correction to the interparticle potential due to a 
finite density of particles may be calculated quite easily 
within a mean field approximation pq . Augmenting the 
free energy F c with an external potential V ex t(r) which 
couples to the concentration, the equilibrium concentra- 
tion is obtained by demanding SF C /Sc = 0. This gives 



c(r) = cq exp 



-V ext (r) 



J d 3 r'c(r')V(r - r') 



(8) 



The mean (or effective) potential is defined by c(r) 



c e" 



-U(r) 



Using this to eliminate the density from the 



previous equation and identifying the external potential 
with that of a fixed particle at the origin, we obtain the 
following self-consistent equation for U(r): 



U(r) 



V(r) + c J d 3 r'V(r - r') exp [- 



£/(r')] • (9) 



Assuming that U <C 1 we can linearise the above inte- 
gral equation. Neglecting the term corresponding to the 
background concentration and Fourier transforming we 
find 



1 2 



1 + c V(q) (B 



C - c Q a 2 )ql 



(10) 



The mean potential is of the same functional form as the 
"bare" interparticle potential, but with a "renormalised" 
characteristic length X 3 (a) — K 3 /(B + C — c a 2 ). Re- 
calling that a similar calculation with a Coulomb poten- 
tial V(q) ~ q~ 2 leads to charge screening pq |, we must 
conclude that, despite the long-ranged interaction, many- 
particle effects arc quite modest in our case. 



B. Correlation Functions 

In the absence of interactions the particles in our 
model have only configurational entropy and their two- 
point correlations are like those of an ideal gas: < 
5c(r)6c(r') >^ S(r — r'). However, elastic interactions 
lead to quite complex correlations, as we show below. 
Defining the dimensionless structure factor 5(q) as the 
Fourier transform of < 6c(r)8c(r') > we find from the 
effective free energy in (|) that 



S(q) 



1 



1 + c V(q) 



= 1 



1 2 

Co a ' ql 



(B + C-c a 2 )ql+K 3 qi 



(11) 



The Fourier transform of the pair correlation function 
h(q) = S(q) — 1 is proportional to —U(q), and therefore 
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the pair correlations themselves are well described by the 
contour plot of Fig.[l] with all the signs reversed. The cor- 
relations are negative along the z and _L axes, indicating 
a reduction in the particle concentration compared to 
the mean value of cq, while in the region bounded by 
the paraboloids z 2 ~ \3Xj_ the correlations are positive, 
indicating an enhancement of the concentration. This 
is broadly consistent with our earlier picture of particle 
chains oblique to the columns. It may be possible to 
measure the particle structure factor by small-angle X- 
ray scattering, and the form of S(q) above would then 
provide a test of the particle-column coupling assumed 
in our theory. 

While the column displacements induce particle corre- 
lations, the particle aggregates in turn modify the col- 
umn fluctuations. In particle-doped lamellar phases, this 
feedback mechanism has been shown to change the usual 
Caille ]37j form of the structure factor pl[ . To calcu- 
late the column displacement correlations in the pres- 
ence of particles, we first define the displacement field 
such that (ViUj) = for a mean concentration Co- This 
requires that we modify the particle-column coupling to 
6c(r)E(u). Next, we define an effective free energy F u for 
the column fluctuations (in a way similar to the definition 
of F c ) by integrating out the concentration field. Carry- 
ing out the integration by steepest descent we obtain the 
effective free energy 



cq (aVi • u e 



-qVj 



')] ■ (I 2 ) 



This shows that the coefficients of all symmetry-allowed 
terms which contain Vj_ • u are "renormalised" by the 
presence of particles. Within our harmonic description 
the only relevant term, obtained by expanding the ex- 
ponential term to quadratic order, is a 2 (V± ■ u) 2 . Re- 
taining only this term in ( |l2[ ) and defining the longitudi- 
nal and transverse parts of u q as q 2 u^ = q(q ■ u q ) and 
ut = u q Ui we obtain the Ulct) ~ U L(T) correlations 
as 



Cq 2 ± + K 3 qi 



(13) 



(B + C- c a 2 )q 2 ± + K 3 qt 



(14) 



The transverse correlations are identical to those of the 
pure columnar phase, which indicates that there are no 
corrections to C and K 3 due to the particles. The longi- 
tudinal correlations while identical in form, have a modi- 
fied bulk modulus B c = B — c a o?. Particles thus produce 
a softening of the columnar phase. 




FIG. 2: A "tail" defect configuration of topological charge 
m — +1 with the thick lines showing the column deformations 
column near the defect. The particle-defect potential vanishes 
on the paraboloids z 2 ~ X3XX (dashed lines), beyond which 
it is mainly repulsive for z > and attractive for z < 0. The 
thin lines show typical trajectories of overdamped particles: 
they are expelled from the regions of dilation (2 > 0) and 
attracted to the regions of compression (z < 0). 



III. PARTICLE-DEFECT INTERACTIONS 

A. Point defects in the columnar phase 

Topological defects are common in liquid crystalline 
phases, being easily produced under conditions of pro- 
cessing such as shear flows. Much of the interesting macro- 
scopic behaviour of liquid crystalline colloids comes from 
the interaction of the colloidal_particles with the topo- 
logical defects of the host 



sum 



For example, the 

solid-like response of particle-doped cholesteric || and 
smectic [Q phases is due to the formation of a stable net- 
work of defects. The colloidal particles play the crucial 
role of stabilising this network, which otherwise would 
anneal away. It is likely that elastic interactions do not 
play a dominant role in such processes, and that short- 
range forces, not captured in a coarse-grained description 
such as ours, are more important. Nevertheless, elastic 
interactions between particles and topological defects are 
known to have important consequences |36] , and so in 
what follows, we calculate this interaction for a point de- 
fect configuration in the columnar phase. The statistical 
mechanics of thermally activated point defects [^7], |3[| 
may be modified in the presence of particles and a knowl- 
edge of particle-defect interactions could be potentially 
useful in a study of this phenomenon. An interesting 
duality (see below) between the columnar and lamellar 
phases has been noted by de Gennes and Pros t P3fl and 
explored in detail by Selinger and Bruinsma [|28| . This 
suggests that it may be of interest to compare particle- 
defect interactions in the lamellar and columnar phases 
and we conclude this section with such a comparison. 

The defect configuration we consider is a single col- 
umn ending p8| , ^fl , illustrated in Fig. |^. This is not 
unduly restrictive, since any linear defect which consists 
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of column endings (eg. a transverse edge dislocation) 
can be constructed out of a collection of such point de- 
fects p3|. The ease with which a column can break de- 
pends on what constitutes the "column", and therefore, 
on what the particular realisation of the columnar phase 
is. In polymeric columnar phases, the columns consist 
of long chains arranged end-to-end along their length. 
Here, thermally activated column scission can occur only 
at those places where one end of a chain meets that of 
another. (Polymeric liquid crystals are found predom- 
inantly in the nematic and isotropic phases; however, 
longer and stiffer biopolymers like DNA can also posi- 
tionally order into columnar phases Since columns 

in discotic and micellar phases are composed of weakly 
bonded molecules, thermally activated scission can take 
place anywhere along the length of the columns. In all 
of these examples, the column scission energy is finite 
and thus a state of thermal equilibrium at non-zero tem- 
perature will always contains a finite number of such de- 
fects. The ends of a scissioned column may either sep- 
arate along the z-axis to form a string of vacancies, or 
may translate past each other in the _L direction to form 
a string of interstitials |2?j . Such strings will begin from 
and end on the column extremities, the lower of which, 
by convention, is called a "tail" (Fig. |J) and the upper a 
"head" |27j . At high temperatures the entropy of wan- 
dering of the string can overcome the energy due to its 
line tension, leading to a string proliferation, and pro- 
ducing isolated "heads" or "tails" |^9|, |3C| ]. Signatures 
of such defects have been observed in recent experiments 
p2[ . An isolated semi-infinite string running along the 
z-axis with its terminus at the origin imposes the con- 
straint that |2l 



Po 



(15) 



where the topological "charge" m is +1 (—1) for a tail 
(head), and po is the areal density of the columns. Clearly 
this is nothing but a line of centers of dilation and is 
easily treated in the framework of the previous section. 
Formally, a column terminus at the origin is identical to a 
line of point particles distributed with the concentration 



(16) 



Combining this with the expression in (§) for the dis- 
placement produced by a distribution of particles, the 
displacement due to a column ending is obtained as J33| 



u 



iqj 



+ A§g| p iq z 



(17) 



The displacement for a collection of column endings dis- 
tributed with the topological "charge" density m(r) fol- 
lows from the superposition of the individual displace- 
ments. Stated mathematically, this means that the total 
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FIG. 3: Plot of the equipotentials of V am (r) ~ V*G(r). The 
coordinates are in units of A3. The contour lines have been 
drawn on equal logarithmically spaced intervals from .001 to 
100, and -0.001 to -100. The potential is odd in z. 



displacement is 



(18) 



with m q being the Fourier transform of m(r). To cal- 
culate the particle-defect interaction, the displacement 
field is decomposed as u q = u q (c q ) + u q (m q ), where 
the first part (due to the particles) is obtained from 
@) and the second part (due to the defects) is ob- 
tained from (|l8|). The free energy then separates out 



as F 



F„ 



F„ 



F cm , with 



and 



p 



2 7 q phlql + XUXl 



iaq\ 



Poq z q'i + \\q A ' 



(19) 



(20) 



The Fourier transforms of the defect-defect and particle- 
defect interactions can now be read off as the kernels of 
TO q m_ q and c q m_ q in ( ^9| ) and (|2Cj ) respectively. The 
transformation to real space can be simplified by noting 
that 



(21) 



Using this property, q\ can be replaced by A§gf in the 
numerators of the kernels in ( fl9| ) and (|2^). Alternatively, 
a conjugate Green's function can be introduced as in |2q| 
and j35[ . With either of the methods, the final form of 
the interaction kernels turn out to be 



V mm {q) = 



■2 2 , \2 4 • 



(22) 
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and 



V cm {q) 



Pa q{ + \iq* 



(23) 



Apart from the line tension terms which we ignore, the 
defect-defect interaction is identical to the result of |H] 
and has been discussed in detail there |34| . The particle- 
defect interaction can be represented in real space as 



V cm (v) = 3 -V 3 z G(x x ,z). 

Pa 



(24) 



As might have been anticipated, the particle-defect po- 
tential is non-central and radially non-monotonic, and is 
odd in the z-coordinate. The behaviour in the z and _L 
directions is simple. Along the columns (for z 2 3> iX3X±) 
the potential decays as sgn(z:)|z| _4 indicating a repulsive 
force for z > and an attractive force for z < 0. In the 
x-y plane, the _L component of the force vanishes by sym- 
metry and the z component goes to zero as Xj°^ 2 . The 
detailed form of the interaction for positive z is shown 
in the equipotential plot of FigJ|. The general conclu- 
sion is that independent of the sign of the defect, the 
particles are attracted to the regions of compression, and 
repelled away from the regions of dilation that the de- 
fect alone would have produced, see Fig.^j. Driven by 
the attractive component of the interaction, and aided 
by Brownian motion, a cloud of particles should collect 
around a column ending, a situation reminiscent of the 
Cottrell clouds of vacancies and interstitials around dis- 
locations in metals |56|. Normally, such clouds decrease 
the dislocation mobility |36| and the same can be as- 
sumed to happen here. Since defect motion is one of the 
main mechanisms of stress relaxation in ordered phases, 
less mobile defect strings may lead to an increase in the 
viscosity of the columnar phase. 

It is important to note that expression can be used 
to obtain the interaction of a particle with any linear 
defect which can be represented as a collection of column 
endings. For example, the interaction of a transverse 
edge dislocation with a particle can be obtained from 
( |l8| ) by setting m(r) ~ S(y)S(z) and repeating the steps 
leading to (|2C 



B. Line defects in the lamellar phase 

There is an interesting line-plane duality p3| , ^| be- 
tween defect-free columnar and lamellar phases which is 
apparent in a pictorial representation: a set of lines repre- 
senting the columns (aligned, say, along the z-axis) when 
rotated by tt/2 can equally well be viewed as a set of 
planes representing the lamellar layers (now stacked nor- 
mal to the z-axis). Mechanically, the lamellar phase is 
a 2d liquid and a Id solid, while the columnar phase is 
a Id liquid and a 2d solid. Though this duality is lost 
in the presence of an arbitrary ensemble of defects (for 



example screw dislocations) , nonetheless there is one im- 
portant exception. It is easy to see from Fig. || that 
an endpoint defect (the terminus of a line) in the colum- 
nar phase corresponds to an edge dislocation (the ter- 
minus of a plane) in the lamellar phase. Going beyond 
pictures, there are close connections in the energetics of 
these "dual" defects, which have been explored in p8[ . 
This prompts us to ask if similar connections may be 
found in the interactions of endpoints and edge disloca- 
tions with particles. To answer this question, we first 
examine the particle-dislocation interaction for an arbi- 
trary distribution of dislocations, and then specialise to 
the case of an edge dislocation. 

The free energy density of a particle-doped lamellar 
phase where the particles couple to the local layer dilation 
has been obtained in Jl9] |. Excluding the entropic terms, 
in our notation it is 

B f 

F = — / d 3 r [{V z uf + A 2 (Viu) 2 + ac(r)V,w] , (25) 

where u is the local displacement of the layers from their 
equilibrium position, B is the bulk modulus, K = B\ 2 
is a second order elastic constant corresponding to layer 
curvature, and the positive constant a is the strength 
of the particle-layer coupling. In the presence of a dis- 
location density b(r) the distortion v = Vu acquires a 
singular transverse part determined by |is| 



V x v = b(r). 

The total distortion in Fourier space is given by 

Vq(b q ) = + tqhq, 



(26) 



(27) 



where /i q , the longitudinal part of the distortion, is ob- 
tained from the equation of equilibrium 5F/6u — 0. The 
solution 



aiq z Cq q z (l - A 2 ^) (z • q x b q ) 



X 2 q\ q 2 + X 2 qi 



q 



(28) 



is a sum of the particle and defect contributions. When 
the distortion in (|27|) is inserted into the (Fourier trans- 
formed) free energy (|2§|), the latter separates into three 
pieces: the first contains only the dislocation density [ p6[ , 
the second only the concentrations [jl9) , and the third is 
a cross-term which represents the particle-dislocation in- 
teraction. Displaying only this last term, we have 



B 



a\ 2 gl 
q 2 + X 2 q\ 



<(z • q x b q )c_ c 



(29) 



This gives the interaction energy of an arbitrary distri- 
bution of particles with an arbitrary number of disloca- 
tions. Two important special cases are linear edge and 
screw dislocation of infinite extent. The dislocation den- 
sity for a screw dislocation of infinite extent and Burg- 
ers vector bz is b(r) = bS(x)8(y)z, which implies that 
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FIG. 4: An infinite edge dislocation of Burgers charge b = +1 
running along the y-axis with the thick lines showing the 
layer deformations near the defect. The dashed lines show 
the parabolas x 2 = 4A|z| for arbitrary values of y. The thin 
lines show typical trajectories of overdamped particles: they 
are expelled from the regions of dilation (x < 0) and attracted 
to the regions of compression (x > 0). Notice that this con- 
figuration is obtained by rotating Fig.[] through n/2. 



b q = 27rb5(q z )z. Inserting this into fl29| ) we find that 
there is no elastic interaction between a particle and the 
screw component of a dislocation within a linear elastic 
theory. For the case of interest, an infinite edge dislo- 
cation running along the y-axis with Burgers vector bz, 
b(r) = bS(x)S(z)y, and b q = 2nbS(q z )z, and so together 
with (|2|) the particle-edge dislocation interaction kernel 
is obtained as 



V cb (q) 



9± 



Transforming to real space, the interaction potential is 



2irbiq x S(q y ). 



(30) 



V cb (v) 



exp 



V^W V4AM 



(31) 



The interaction is long-ranged along the z-axis (for 
x 2 <C 4A|z|) and exponentially damped along x-axis (for 
x 2 4A|z|). By symmetry, the z-component of the force 
vanishes on the z-axis, while both the x and z com- 
ponents are exponentially small near the £-axis. Irre- 
spective of the sign of 6, the particles are attracted to 
the compression regions that the dislocation alone would 
have produced, and repelled away from the dilation re- 
gions, see Fig.([|). To compare this interaction with its 
dual in the columnar phase, we first note that the in- 
teraction is characterised by the parabolas x 2 — 4A|z| 
which occur frequently in problems related to disloca- 
tions in lamellar phases f23fl . The particle-defect poten- 
tial in the columnar phase is characterised instead by 
the paraboloids z 2 ~ \-$x±. In both cases the maxi- 
mum elastic deformation is confined to the regions within 
these surfaces, compare Fig.(^) and Fig.(Q). The force is 
long-ranged along the solid-like directions (the _L direc- 
tion in the columnar phase, and the z direction in the 



lamellar phase) and short-ranged in the liquid-like direc- 
tions. There is no interaction between a particle and an 
edge dislocation when both have the same z coordinate, 
a feature also seen in the interaction between two edge 
dislocations (^3). In contrast, the particle-defect interac- 
tion in the columnar phase is strongly repulsive in the 
z direction. The formation of Cottrell clouds, and the 
consequent reduction in the dislocation mobility should 
be expected in the lamellar phase as well. 



IV. CONCLUSION 

To summarise, we have calculated the elastic interac- 
tions of colloidal particles in a columnar hexagonal phase 
by modelling them as centers of dilation. The interaction 
is long-ranged and non-central, and has both attractive 
and repulsive parts. We have shown that this interac- 
tion leads to correlations between the particle positions, 
and to a reduction in the bulk modulus B of the colum- 
nar phase. We have also shown how particles interact 
elastically with defects in columnar and lamellar phases. 
In both phases, particles aggregate around the defects 
driven by the attractive component of the interaction and 
reduce the defect mobility. One likely consequence of this 
is an increase in the viscosity of the phase. 

We mention below several of the directions in which 
the present work may be extended. First, the elastic in- 
teractions of particles of more general shapes like cylin- 
ders and discs should be examined. Such inclusions can 
no longer be modelled as force dipoles, and force distri- 
butions containing higher multipoles will be needed [p4[ . 
The interaction of plate-like particles would have a direct 
relevance to the physics of clay-block copolymer compos- 
ites p2[ . Second, the effects of boundaries should be 
investigated. Free boundaries require the introduction 
of image sources, which will have long-range interactions 
with particles in the bulk. Thus, boundary effects can be 
expected to strongly influence all processes in the bulk, 
and in particular, the formation of particle aggregates. 
Third, the nature of the topological defects surrounding 
a particle should be examined. Our coarse-grained de- 
scription cannot predict the defect configuration (though 
useful information can be gained if we assume it) and a 
more microscopic description would be better suited for 
the purpose. 

As a final, somewhat speculative remark, we should 
indicate that the present calculation may be of some 
relevance in biological systems where inclusions may be 
present in tissues of columnar hexagonal structure [ 38 1 . 
Biological inclusions generate local inward forces |ll], 39 1 , 
and a minimal description of such inclusions could be in 
terms of the centers of compression (corresponding to a 
negative a) in our theory. Similar elastic interactions be- 
tween proteins in a lamellar phase p0| , and cells in an 
isotropically elastic tissue have already been investigated 
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